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Recent  progress  in  thermoelectric  power  production  using  Bismuth  Telluride  Bi2Te3  semiconductor 
modules  has  revealed  the  potential  to  effectively  convert  large  volumes  of  low  temperature  industrial 
waste-heat  to  electricity.  In  order  to  render  the  process  more  cost  effective,  greater  understanding  of 
the  effects  of  external  influences  on  the  module’s  power  output  is  necessary.  Such  an  understanding 
would  facilitate  the  design  of  thermoelectric  generators  which  serve  to  exploit  available  waste-heat. 
To  this  end,  an  experimental  study  is  performed  on  the  most  adjustable  operating  parameter  on  a  ther¬ 
moelectric  liquid-to-liquid  generator,  the  electrical  load  resistance.  A  test  stand  apparatus  is  built  apply¬ 
ing  a  temperature  gradient  on  commercially  available  Bi2Te3  thermoelectric  modules  by  means  of  an 
injection  and  a  rejection  of  heat  brought  upon  by  counter  current  hot  and  cold  liquids.  The  thermoelectric 
power  production  relative  to  an  increasing  electrical  load  is  investigated  by  means  of  an  analysis  of 
experimentally  measured  results  in  which  the  thermal  input  conditions  are  varied.  The  results  detail 
the  thermoelectric  characteristics  of  a  liquid-to-liquid  generator  under  an  increasing  electrical  load  resis¬ 
tance  by  identifying  the  optimal  electrical  load  resistance  for  peak  thermoelectric  production.  A  correla¬ 
tion  between  peak  thermoelectric  power  and  thermal  input  conditions  is  presented  as  well  as  an 
investigation  into  the  validity  of  electrical  load  matching. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  modules  are  now  being  considered  in  a  variety 
of  power  generation  applications  in  which  the  caloric  resource  is 
waste-heat  from  processes  such  as  pulp  and  paper  production,  data 
centre  cooling  systems  and  engine  exhaust.  In  such  considerations, 
the  goal  is  to  generate  electricity  from  low-grade  waste-heat  that  is 
within  a  temperature  range  of  25-225  °C.  In  this  range,  the  semi¬ 
conductor  Bismuth  Telluride  (Bi2Te3)  has  been  shown  [1-3]  to  be 
the  most  effective  commercially  available  thermoelectric  material 
for  exploiting  the  thermoelectric  phenomenon  known  as  the  See- 
beck  effect.  For  this  reason,  many  studies  on  low  temperature 
applications  of  thermoelectric  power  generation  investigate  the 
properties  of  Bi2Te3  thermoelectric  modules  [4-6]. 

Generally,  thermoelectric  power  generation  aims  to  exploit  low 
cost  heat  sources  [7-9].  For  example,  [10-12]  detail  the  potential 
use  of  thermoelectric  modules  in  automotive  waste-heat  recovery, 
[13]  demonstrate  the  use  of  thermoelectric  materials  in  converting 
biomass  cookstove  waste-heat  to  electricity,  and  [14]  investigate 
the  harnessing  of  solar  energy  by  coupling  thermoelectric  modules 
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with  concentrator  lenses.  In  such  applications,  variable  operating 
conditions  imply  variable  thermal  input  conditions.  It  is  commonly 
argued  that  such  dynamic  systems  require  a  varying  electrical  load 
for  optimal  power  output;  indeed,  [15]  describe  the  importance  of 
electrical  load  matching  under  dynamic  thermal  conditions  in 
exhaust  waste-heat  recovery  applications.  Similarly,  [16]  discuss 
the  importance  of  electrical  load  matching  due  to  changing  ther¬ 
mal  conditions  for  a  generic  low  temperature  thermoelectric 
power  generation  system.  This  entails  that  an  application  in  which 
the  working  conditions  vary  greatly,  such  as  an  automotive 
exhaust  recovery  application  or  a  solar  energy  harvesting  appara¬ 
tus,  requires  a  varying  electrical  load  resistance  shifting  according 
to  fluctuations  in  thermal  conditions  in  an  effort  to  maximize  the 
thermoelectric  power  output.  The  adjustment  in  the  electrical  load 
resistance  that  is  necessary  when  thermal  conditions  change  is 
attributed  to  electrical  load  matching  with  the  internal  electrical 
resistance  of  the  generator  and  is  referred  to  as  Maximum  Power 
Point  Tracking  (MPPT).  Indeed,  by  using  Bi2Te3  modules  to  convert 
low  grade  thermal  energy  from  liquid  metal  to  electricity  and  by 
performing  the  same  test  under  two  different  loads,  [17]  illustrated 
that  the  electrical  load  influences  the  output  voltage  in  thermo¬ 
electric  power  generation.  This  implies  that  there  exists  an  optimal 
load  resistance.  However,  under  varying  thermal  conditions,  [17] 
did  not  adjust  the  electrical  load  for  peak  power  production  and 
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Nomenclature 

Symbol 

Description  (Unit) 

opt 

optimal  (-) 

A 

contact  area  of  a  thermocouple  pellet  (m3) 

P 

power  (W) 

H 

hot  side  (-) 

Q 

heat  transfer  rate  (W) 

C 

cold  side  (-) 

Rl 

electrical  load  resistance  (Q) 

cP 

specific  heat  at  constant  pressure  (J/kg  I<) 

Rint 

internal  electrical  resistance  (Q) 

I 

electrical  current  (A) 

T 

temperature  (K) 

k 

thermal  conductivity  (W/m  K) 

V 

voltage  (A) 

L 

length  of  a  thermocouple  (m) 

V 

fluid  velocity  (m/s) 

max 

maximum  (-) 

$p,n 

Seebeck  coefficient  (V/K) 

N 

number  of  thermocouples  (-) 

P 

electrical  resistivity  (Q  m) 

oc 

open  circuit  (-) 

therefore  did  not  examine  the  influence  of  a  variable  electrical 
load. 

The  purpose  of  this  study  is  to  contribute  to  the  development  of 
MPPT  systems  by  evaluating  the  power  versus  electrical  load  re¬ 
sponse  of  a  novel  liquid-to-liquid  thermoelectric  generator  with 
respect  to  its  thermal-hydraulic-electrical  performance.  To  this 
end,  a  reliable  test  facility  for  an  experimental  characterization  of 
the  effects  of  an  increasing  load  resistance  on  thermoelectric 
power  output  is  established.  A  temperature  gradient  across  a  series 
of  forty  thermoelectric  modules  is  produced  by  the  flow  of  hot  and 
cold  water.  In  particular,  this  study  shows  that  under  an  increasing 
load  resistance,  changing  thermal  input  conditions  have  negligible 
effects  on  the  optimal  electrical  load  resistance,  noted  flLi0pt,  for 
power  output.  Secondly,  a  correlation  is  developed  relating  the 
generator’s  peak  thermoelectric  power  output  to  thermal  input 
conditions;  the  correlation  shows  that  peak  power  displays  expo¬ 
nential  behavior  relative  to  temperature  and  root  function  behav¬ 
ior  relative  to  flow  rates.  Thirdly,  it  is  shown  that  thermoelectric 
power  is  less  sensitive  to  electrical  load  error  when  the  load  is 
greater  than  optimal  rather  than  less  than  optimal.  Finally,  the 
present  work  demonstrates  that  the  RLt0pt  for  thermoelectric  pro¬ 
duction  under  an  increasing  electrical  load  resistance  is  strictly  less 
than  the  internal  resistance  of  the  generator  for  all  thermal  input 
conditions  tested. 

2.  Thermoelectric  power 

The  focus  of  this  study  is  the  thermoelectric  power  production 
of  a  generator  relative  to  an  increasing  electrical  load  resistance. 
The  relationship  between  thermoelectric  power  and  its  electrical 
load  is  detailed  in  [18,19]  in  which  a  thermoelectric  module  made 
up  of  N  thermocouples,  such  as  in  Fig.  1,  is  considered. 

Assuming  constant  and  isotropic  material  properties  and  an 
interconnecting  material  Seebeck  coefficient  of  zero,  the  one¬ 
dimensional  heat  equation  with  boundary  conditions  T(x  =  0 )  =  TC 
and  T(x  =  L)  =  TH  yields  the  temperature  profile: 


T 


12  p*2 


2kAz 


Te¬ 


rn 


in  which  the  hot  and  cold  side  temperatures  are  noted  TH  and  Tc 
respectively,  A T=TH  -  Tc,  and  I  is  the  electric  currant.  Energy  con¬ 
servation  at  the  hot  junction,  x  =  L,  requires  that  the  heat  flux  “in” 
be  equal  to  the  heat  flux  “out”  plus  the  energy  absorbed  by  the  See¬ 
beck  effect: 


Qh  —  Qout  ~  IXp,nTH.  (2) 

In  Eq.  (2),  ap,n  =  otp  -  an  is  the  difference  in  Seebeck  coefficients 
of  the  p- type  and  n-type  semiconductor  pellets. 

Combining  Eqs.  (1)  and  (2)  with  Fourier’s  law  of  conduction 
qout  =  k(2A)  g  at  x  =  L,  yields  the  following  energy  balance: 

2kA  l2  pi 

qh  =  -j-Ar--£-  +  /ap,nrH.  (3) 

Similarly,  at  x  =  0,  the  cold  side  heat  flux  is: 

2 kA  l2pl  .  _  ... 

Qc  —  AT  +  — +  IoCpnTc .  (4) 

The  power  of  a  single  module  is  the  sum  of  the  heat  flux  differ¬ 
ences  of  the  N  thermocouples: 

P  =  N(qh  -  qc)  =  -I2  +  Map,„AT.  (5) 

Since  the  term  2 NpL/A  represents  the  sum  of  the  electrical 
resistances  of  the  thermocouples,  noted  the  internal  resistance  Rint, 
Eq.  (5)  can  be  expressed  as: 

I2RL  =  -I2Rin[  +  NIappAT.  (6) 

Solving  for  a  nonzero  electrical  current  in  Eq.  (6)  and  combining 
with  Eq.  (5)  yields  the  power  output  in  terms  of  the  temperature 
difference  and  the  electrical  load  resistance: 


p  =  (NATaMf  (NATxP:n)2RiM 
Rint  +  RL  (Rj„t  +  Rl)2 


(7) 


An  analysis  of  the  above  equation  easily  shows  that  the  maxi¬ 
mum  power  output  occurs  when  the  electrical  load  resistance  im¬ 
posed  on  the  system  is  equal  to  the  internal  electrical  resistance  of 
the  generator;  that  is  to  say,  peak  power  is  obtained  when  the  elec¬ 
tric  load  normalized  by  the  internal  resistance  attains  unity: 

Rl  =  RL/Rint  =  l-  (8) 


This  result  yielding  a  peak  power  of: 

D  _  (NATaPi„)2 

4  Rl 


(9) 


Fig.  1.  Thermoelectric  production  from  thermocouples  electrically  in  series  and  *s  often  used  when  modeling  peak  thermoelectric  power  Such  as  in 
thermally  in  parallel.  the  dynamic  thermal  conditions  studied  in  Eakburanawat  and 
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Boonyaroonate  [20]  and  Crane  and  Bell  [15].  The  validity  of  load 
matching  for  peak  thermoelectric  power  is  investigated  in  the  con¬ 
text  of  an  increasing  electrical  load. 


3.  Experimental  setup 

An  experimental  test  stand  is  built  in  order  to  measure  power 
output  of  thermoelectric  modules  with  an  increasing  electrical 
load  resistance.  The  test  stand  supplies  heat  by  channelling  hot 
water  between  two  layers  of  thermoelectric  modules  that  are  elec¬ 
trically  in  series  while  drawing  heat  away  from  the  outer  side  of 
the  module  layers  with  cold  water. 

The  central  piece  of  the  test  stand  is  a  generator  consisting  of 
two  sets  of  twenty  commercially  available  Bi2Te3  thermoelectric 
modules  TEG2-07025HT-SS  placed  on  either  side  of  a  hot  fluid  flow 
channel.  Two  cold  fluid  flow  channels  enclose  the  system  dissipat¬ 
ing  heat  thereby  stabilising  a  temperature  gradient  across  the 
modules.  The  result  is  an  aluminum  encased  purpose  built  thermo¬ 
electric  liquid-to-liquid  generator  assembled  by  Thermal  Electron¬ 
ics  Corp.,  and  presented  in  Fig.  2. 

In  this  way,  the  contrast  in  the  hot  temperature  field  with  the 
cold  temperature  field  generates  a  thermoelectric  current  known 
as  the  Seebeck  effect.  A  thermoelectric  module  is  made  up  of  an 
assemblage  of  p-n  junctions  that  are  electrically  in  series  due  to 
the  alternating  n-type  and  p- type  Bi2Te3  semiconductor  pellets. 
The  module  assemblage  is  held  in  place  with  ceramics,  is  thermally 
in  parallel  and  measures  40  mm  x  40  mm  x  3.75  mm. 

The  generator  in  which  the  modules  are  encased  measures 
47.22  cm  by  9.07  cm  by  1.73  cm.  The  hot  and  cold  fluid  flow  tubes 
have  an  interior  diameter  of  9.52  mm  and  an  exterior  diameter  of 
12.70  mm  and  water  is  used  as  the  working  fluid.  The  hot  water 
pipes  of  the  generator  are  connected  to  a  MicroTherm  CMX  Series 
Temperature  Control  System  which  is  a  combined  pump  and  heating 
closed  loop.  Thermocouple  temperature  sensors  are  positioned  at 
each  fluid  entry  and  exit.  Pressure  gauges  are  placed  at  the  entry 
and  exit  of  the  hot  fluid  flow.  Data  acquisitioning  of  fluid  temper¬ 
ature,  fluid  pressure,  electrical  load  resistance  and  voltage  are  per¬ 
formed  with  DataStudio  software  in  which  sensor  readings  are 
relayed  to  the  computer  via  the  interface  ScienceWorkshop  750. 
A  carbon  wire  rheostat  with  a  manual  cursor  is  used  to  vary  the 
electrical  load  resistance  from  0  to  40  Q.  Flow  turbulating  inserts 
are  placed  in  all  fluid  flow  channels  in  an  effort  to  improve  thermal 
transport  between  the  thermoelectric  materials  and  the  water.  In¬ 
deed,  [21  ]  identify  in  a  previous  study  that  the  flow  turbulating  in¬ 
serts  with  alternating  tabs  every  7.9  mm  on  a  9  mm  x  0.5  mm 


Fig.  2.  Aluminum  encased  thermoelectric  liquid-to-liquid  generator. 


shim  greatly  improve  power  output  over  a  large  range  of  thermal 
conditions.  The  control  system  of  the  experimental  apparatus  is 
presented  in  Fig.  3.  It  is  important  to  note  that  the  hot  and  cold 
fluid  flow  channels  are  in  opposing  directions.  This  is  noteworthy 
since  it  is  also  common  -  in  an  effort  to  increase  the  quantity  of 
thermoelectric  material  used  -  to  stack  the  modules  in  many  layers 
as  was  done  by  [22].  However,  this  leads  to  alternating  parallel  and 
countering  hot  and  cold  fluid  flow  channels  which  differs  mechan¬ 
ically  from  the  test  stand  of  the  present  work.  In  the  present  study, 
parallel  flow  is  purposely  avoided  in  order  to  favor  a  homogeneous 
temperature  gradient  the  length  of  the  generator. 

The  electrical  resistances  RB ,  Rlf  and  R2  identified  in  Fig.  3  are 
measured  to  be  1.0  Q,  39.0 1<Q,  and  14.8  1<Q  respectively  within 
an  accuracy  of  0.8%  using  a  BI<  Precision  Digital  Multimeter  Model 
388-B.  The  voltage  across  the  TEG  is  measured  using  a  voltage  di¬ 
vider  setup  such  that: 

VTeg  =  VaR^A.  (10) 

The  thermoelectric  power,  noted  P,  is  calculated  by  measuring 
VA  and  VB  and  coupling  them  with  the  measured  electrical  resis¬ 
tances  of  the  circuit  such  that: 


P  =  VtecI  =  Vtectt. 

Kb 


(ii) 


Similarly,  the  electrical  load  resistance  is: 

r,  VtEG  t  /  Rb 

Kl  =  — r  =  V  TEG  YT- 

1  v  B 


(12) 


The  voltages  VA  and  VB  are  measured  using  Science  Workshop 
Interface  to  an  accuracy  of  ±10  mV  and  ±3  mV  respectively.  In  this 
study,  the  thermoelectric  power  output  readings  and  the  electrical 
load  resistance  readings  are  calculated  to  each  have  a  maximum 
uncertainty  of  3.7%  for  all  thermal  input  conditions  tested. 


4.  Results  and  discussion 

The  following  test  cases  characterize  the  effects  of  an  increasing 
electrical  load  resistance  on  the  thermoelectric  production  of  a  li¬ 
quid-to-liquid  generator  for  varying  thermal  input  conditions.  Par¬ 
ticular  attention  is  placed  on  the  electrical  load  resistance  that 
yields  peak  power,  noted  RL,opt-  As  such,  the  present  study  focuses 
on  the  influence  that  thermal  conditions  have  on  the  RLt0pt  under  a 
variable  electrical  load  resistance. 

In  order  to  identify  peak  power  with  respect  to  the  electrical 
load,  the  power  output  is  normalized  by  the  maximum  power  out¬ 
put  of  the  specified  test  case  such  that: 

P*  —  P/P  max  •  (13) 

In  this  way,  the  maximum  point  on  the  normalized  power  out¬ 
put  curves  for  each  test  case  attains  unity  and  optimal  electrical 
loads  are  easily  compared. 

In  order  to  investigate  the  thermoelectric  effect  under  a  variable 
electrical  load,  power  output  of  the  forty  module  liquid-to-liquid 
thermoelectric  generator  is  measured  under  different  thermal  in¬ 
put  conditions.  These  conditions  are  varied  in  three  ways:  firstly, 
by  changing  the  temperature  difference  between  the  hot  and  cold 
fluids;  secondly,  by  changing  the  flow  rate  of  the  hot  fluid;  and 
thirdly,  by  changing  the  flow  rate  of  the  cold  fluid. 

4.1.  Thermal  input  variation  due  to  temperature 

The  first  set  of  experiments  examines  the  effects  of  a  change  in 
thermal  conditions  by  varying  the  hot  inlet  temperature.  The  flow 
rates  for  the  hot  and  cold  sides  are  maintained  at  4.37  L/min  and 
5.00  L/min  respectively,  the  cold  fluid  inlet  temperature  remains 
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Table  1 

Inlet  and  outlet  temperatures  for  test  cases  of  Section  4.1. 


ATavg 

8.6  °C 

16.4  °C 

25.3  °C 

33.0  °C 

50.3  °C 

Thm  (°C) 

30.6 

40.6 

48.9 

59.4 

79.4 

Ti.OUt  (°C) 

29.4 

40.0 

46.7 

58.3 

75.0 

TC,in  (°C) 

21.1 

22.8 

21.1 

23.9 

23.3 

Tc.out  (°C) 

21.7 

25.0 

23.9 

27.8 

30.6 

Fig.  4.  Effect  of  an  increasing  electrical  load  resistance  on  power  output. 


Fig.  5.  Effect  of  an  increasing  electrical  load  resistance  on  normalized  power  output 
relative  to  AT. 


between  21.1  and  23.9  °C,  and  the  hot  fluid  inlet  temperature  is  in¬ 
creased  from  30.6  to  79.4  °C.  The  electrical  load  resistance  is  in¬ 
creased  at  a  rate  of  6.9  ±2.1  Q/s.  Each  test  case  generates  a 
temperature  difference  between  the  heat  source  and  the  heat  sink 
and  is  identified  by  the  difference  in  the  average  cold  side  temper¬ 
ature  and  the  average  hot  side  temperature  such  that: 

ATavg  —  (T H,in  +  T/J,out)/2  ~  C^C,in  +  Tc,out) /^-  (14) 

It  is  important  to  note  that  the  temperature  difference  across 
the  actual  modules  embedded  in  the  generator  is  strictly  less  than 
A Tavg  due  to  the  thermal  resistance  of  the  forced  convection  and 
the  interface  materials. 

The  temperature  readings  at  the  fluid  flow  inlets  and  outlets  for 
each  of  the  test  cases  of  this  section  are  given  in  Table  1. 

The  results  illustrated  in  Fig.  4  show  that  the  thermoelectric 
power  output  for  all  temperatures  tested  have  similar  trends.  In 
particular,  for  each  test  case,  power  output  increases  with  increas¬ 
ing  electrical  load  resistance  until  a  peak  power  is  reached  at  an 
optimal  load  RL,oPt-  Beyond  ftL>opti  the  thermoelectric  power  de¬ 
creases  more  gradually.  It  is  important  to  note  that  Fig.  4  is  pre¬ 
sented  on  a  Logarithmic  scale  in  order  to  illustrate  the  similar 
trends  amongst  the  different  test  cases  despite  the  large  range  in 
power  outputs. 

In  Fig.  5,  the  power  output  is  normalized  by  the  maximum 
power  obtained  in  each  test  case.  It  is  shown  that  all  temperature 
fields  generate  similar  power  output  trends  relative  to  the  electri¬ 
cal  load.  More  specifically,  the  normalized  power  output  curves  for 
all  temperature  conditions  tested  collapse  into  one  curve  suggest¬ 
ing  that  MPPT  due  to  a  variation  in  thermal  conditions  is  not  nec¬ 
essary  relative  to  the  electrical  load  resistance.  Indeed,  the  power 
output  for  each  thermal  input  condition  is  above  99.0%  when  with¬ 
in  the  range  of  all  RLt0pt. 

In  order  to  illustrate  the  peak  power  and  optimal  load  trends 
relative  to  the  temperature  differential,  the  data  points  of  each  test 
case  are  curve  fit  to  a  seventh  order  polynomial  from  which  the 
coordinates  of  the  maximum  point  are  identified.  The  peak  power 
and  optimal  load  trends  are  illustrated  in  Fig.  6.  Consistent  with  Eq. 
(9),  the  maximum  power  output  increases  with  an  approximate 
parabolic  curve  relative  to  AT.  Conversely,  despite  a  large  change 
in  thermal  input  conditions  and  thermoelectric  power  production, 
the  RLiopt  increases  only  gradually  from  14.8  to  17.2  Cl.  This  is 
attributed  to  an  increase  in  material  internal  resistance  brought 
upon  by  an  increase  in  module  mean  temperature;  such  an  in¬ 
crease  in  internal  electrical  resistance  subsequently  increases  the 
optimal  electrical  load. 
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Fig.  6.  Maximum  thermoelectric  power  output  and  corresponding  optimal  electri¬ 
cal  load  resistance  with  respect  to  AT. 

4.2.  Thermal  input  variation  due  to  hot  flow  rate 

The  effects  of  an  increasing  load  resistance  on  power  production 
for  changing  thermal  input  conditions  in  which  the  hot  flow  rate  is 
varied  are  illustrated  in  Fig.  7.  The  hot  fluid  inlet  temperature  re¬ 
mains  between  78.9  and  81.1  °C  and  the  cold  fluid  inlet  tempera¬ 
ture  remains  between  20.0  and  20.6  °C.  The  cold  fluid  flow  rate  is 
maintained  at  2.5  L/min  for  each  test  case  and  the  hot  fluid  flow 
rate  is  increased  from  0.49  to  4.37  L/min.  It  is  important  to  note 
that  the  cross  sectional  area  of  the  flow  tubes  is  fixed  and  that 
any  variation  in  flow  rate  corresponds  to  a  variation  in  fluid  veloc¬ 
ity.  The  electrical  load  resistance  is  increased  at  a  rate  of 
7.8  ±2.6  O/s. 

The  temperature  readings  at  the  inlets  and  the  outlets  of  the 
flow  fluids  for  the  test  cases  of  this  section  are  given  in  Table  2. 

The  normalized  power  output  curves  of  Fig.  7  show  that  each 
test  case  experiences  a  sharp  increase  in  power  output  until  reach¬ 
ing  a  peak  power  at  an  optimal  electrical  load.  In  particular,  an  in¬ 
crease  in  hot  fluid  flow  rate  is  not  shown  to  have  a  measurable 
effect  on  RL,opt-  More  specifically,  for  the  range  of  RL,opt,  all  test  cases 
are  above  99.6%  power  output. 

Fig.  8  plots  the  peak  power  and  optimal  electrical  load  trends 
for  each  test  case.  The  results  differ  from  the  observations  of  Sec¬ 
tion  4.1  in  that  the  peak  power  curve  exhibits  a  more  linear  or  root 
function  behavior  rather  than  parabolic  while  the  RL,opt  curve 
exhibits  a  constant  behavior.  The  peak  power  increase  is  attributed 
to  an  increase  in  convective  heat  transfer  to  the  embedded  thermo¬ 
electric  modules  thereby  increasing  the  temperature  difference 
across  the  modules.  This  increased  temperature  difference  in- 


Fig.  7.  Effect  of  an  increasing  electrical  load  resistance  on  normalized  power  output 
relative  to  hot  fluid  flow  rate. 


creases  power  exponentially  in  accordance  to  Eq.  (9).  The  peak 
power  linear  to  root  function  behavior  is  attributed  to  a  rate  of  in¬ 
crease  of  thermal  transport  with  respect  to  the  hot  flow  rate  that 
decreases  with  increasing  flow  rates.  For  larger  flow  rates,  the  peak 
power  output  is  expected  to  plateau  due  to  the  fact  that  the  upper 
limit  of  the  hot  side  of  the  modules  is  the  temperature  of  the  hot 
inlet  fluid. 

4.3.  Thermal  input  variation  due  to  cold  flow  rate 

A  third  way  to  vary  the  thermal  input  conditions  of  a  thermo¬ 
electric  liquid-to-liquid  generator  is  to  vary  the  cold  fluid  flow  rate. 
Fig.  9  illustrates  the  effects  of  an  increasing  electrical  load  resis¬ 
tance  on  thermoelectric  power  production  for  different  cold  fluid 
flow  rates.  The  hot  fluid  inlet  temperature  remains  between  78.3 
and  79.4  °C  and  the  cold  fluid  inlet  temperature  remains  between 
19.4  and  20.0  °C.  The  hot  fluid  flow  rate  is  maintained  at  2.13  L/ 
min  for  each  test  case  and  the  cold  fluid  flow  rate  is  increased  from 
2.5  to  10.0  L/min.  The  electrical  load  resistance  is  increased  at  a 
rate  of  7.8  ±  2.6  Q/s.  It  is  important  to  note  that  the  cold  fluid  flow 
rates  are  larger  than  the  hot  fluid  flow  rates  since  the  apparatus  has 
four  cold  fluid  flow  channels  compared  to  two  hot  fluid  flow 
channels. 

The  temperature  readings  at  the  inlets  and  the  outlets  of  the 
flow  fluids  for  the  test  cases  of  this  section  are  given  in  Table  3. 

Fig.  9  illustrates  the  effect  of  an  increasing  electrical  load  on  the 
normalized  power  output  for  these  test  cases.  No  measurable 
change  in  the  RLi0pt  due  to  the  variations  in  cold  inlet  flow  rate  is 
observed. 

The  results  illustrated  in  Fig.  10  shows  that  peak  power  exhibits 
root  function  behavior  with  respect  to  the  cold  fluid  flow  rate.  The 
increase  in  power  output  is  attributed  to  the  increase  in  cold  fluid 
flow  which  improves  the  heat  sink  thermal  transport  effectively 
increasing  the  temperature  difference  across  the  embedded  mod¬ 
ules.  The  maximum  power  output  curve  appears  as  a  root  function 
relative  to  the  cold  fluid  flow  rate  since  the  rate  of  increase  of  ther¬ 
mal  transport  with  respect  to  the  cold  flow  rate  decreases  with 
increasing  flow  rate.  The  results  illustrated  in  Fig.  10  also  show 
that,  despite  an  important  increase  on  the  power  output,  a  change 
in  thermal  conditions  due  to  cold  fluid  flow  rate  decreases  the  RLr 
oPf  This  is  attributed  to  the  increased  cold  fluid  flow  lowering  the 
mean  temperature  of  the  modules  which  decreases  the  internal 
resistance  of  the  materials  thereby  decreasing  the  optimal  electri¬ 
cal  load. 

4.4.  Generalized  trends  and  thermoelectric  peak  power 

The  thermoelectric  power  output  trends  for  all  of  the  thermal 
input  conditions  tested  in  the  previous  sections  are  curve  fit  to 
one  curve  in  order  to  illustrate  its  rate  of  change  relative  to  the 
electrical  load;  the  results  are  illustrated  in  Fig.  11. 

For  all  thermal  input  conditions  tested,  a  sharp  increase  in 
power  output  is  produced  as  the  electrical  load  approaches  RLt0pt 
where  upon  maximum  power  is  attained.  As  the  electrical  load  in¬ 
creases  beyond  RL,opt,  thermoelectric  power  decreases  more  gradu¬ 
ally  at  a  nearly  constant  rate.  This  implies  that  an  MPPT  should  err 
on  the  side  of  an  electrical  load  resistance  that  is  greater  than  RLi0pt 
rather  than  less  than. 

Specifically,  the  larger  jjj£  values  for  RL  <  RLy0pt  imply  that  the 
power  output  is  sensitive  to  small  changes  in  RL  in  this  range.  How¬ 
ever,  for  Rl  >  RLjoPt ,  values  of  ^  near  zero  are  produced  tending  to  a 
constant  value  implying  that  the  thermoelectric  power  output  is 
less  sensitive  to  electrical  loads  that  are  greater  than  optimal. 

In  order  to  factor  in  all  thermal  conditions  into  a  single  non- 
dimensional  term,  x*  is  defined: 
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Table  2 

Inlet  and  outlet  temperatures  for  test  cases  of  Section  4.2. 


Hot  flow  (°C): 

0.49  L/min 

1.02  L/min 

1.60  L/min 

2.13  L/min 

3.40  L/min 

3.86  L/min 

4.37  L/min 

TH,in 

79.4 

80.6 

78.9 

79.4 

78.9 

78.9 

81.1 

TH.out 

57.8 

62.2 

68.3 

71.1 

72.2 

73.9 

75.0 

TC,in 

20.0 

20.0 

20.0 

20.0 

20.6 

20.6 

20.6 

Tc,out 

30.0 

30.0 

32.2 

32.2 

32.8 

32.8 

32.2 

Fig.  8.  Maximum  thermoelectric  power  output  and  corresponding  optimal  electri¬ 
cal  load  resistance  with  respect  to  hot  fluid  flow  rate. 


Fig.  10.  Maximum  thermoelectric  power  output  and  corresponding  optimal 
electrical  load  resistance  with  respect  to  cold  fluid  flow  rate. 


Fig.  9.  Effect  of  an  increasing  electrical  load  resistance  on  normalized  power  output 
relative  to  cold  fluid  flow  rate. 


dP* 


Fig.  11.  Effect  of  an  increasing  electrical  load  resistance  on  normalized  power 
output  and  its  rate  of  change  relative  to  the  electrical  load. 


Table  3 

Inlet  and  outlet  temperatures  for  test  cases  of  Section  4.3. 


Cold  flow  (°C): 

2.5  L/min 

5.0  L/min 

7.5  L/min 

10.0  L/min 

THjn 

79.4 

78.9 

78.3 

79.4 

T H,out 

71.1 

70.0 

69.4 

70.0 

TC,in 

20.0 

20.0 

19.4 

20.0 

Tc.out 

32.2 

25.0 

22.8 

23.9 

X* 


20,200 


Vh 


l  yj CpflAT avg J  l  \J' Cp,C^T  ( 


Vc 


(15) 


in  which  vH  and  vc  are  the  hot  and  cold  side  fluid  velocities. 

The  first  two  terms  on  the  right  side  of  Eq.  (15)  are  Eckert  num¬ 
bers.  Their  numerators  represent  the  kinetic  energy  of  the  fluid 
flows  while  their  denominators  represent  the  hot  and  cold  side  en¬ 
thalpy  differences;  The  last  term  on  the  right  side  of  Eq.  (15)  is  the 
Carnot  efficiency.  It  is  duly  noted  that  the  non-dimensional  term  x* 
of  Eq.  (15)  is  a  root  function  with  respect  to  the  fluid  flow  and  an 
exponential  function  with  respect  to  the  temperature  difference. 


This  non-dimensional  parameter  quantifies  the  thermal  input  into 
the  generator  as  a  function  of  the  hot  and  cold  fluid  temperatures 
and  flow  rates. 

Fig.  12  compiles  the  measured  peak  power  measurements  from 
ten  test  cases.  The  thermal  input  conditions  vary  in  temperature 
difference  between  the  hot  and  cold  sides,  hot  fluid  flow  rate  and 
cold  fluid  flow  rate.  The  peak  power  points  are  presented  as  sets 
of  data  in  which  the  thermal  input  parameters  are  fixed  with 
exception  to  either  the  hot  fluid  flow  rate,  the  cold  fluid  flow  rate 
or  the  temperature  difference.  The  results  show  that  the  power 
output  of  the  tested  generator  remains  within  10%  of  the  value  pre¬ 
dicted  by  the  non-dimensional  term  x*  for  all  thermal  input  condi¬ 
tions  tested.  In  this  way,  the  peak  power  output  of  the  generator 
can  be  predicted  for  given  thermal  input  parameters.  The  results 
agree  with  the  parabolic  behavior  relative  to  temperature  and 
the  root  function  behavior  relative  to  flow  rate. 

Fig.  13  compares  measured  RL<opt  values  to  the  thermal  input 
conditions  in  terms  of  the  non-dimensional  term  x*.  No  measurable 
trends  in  the  RLy0pt  with  respect  to  the  thermal  input  conditions 
tested  are  observed.  In  particular,  the  RL>opt  is  remains  within  15% 
of  15.8  Q. 
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Fig.  14.  Effect  of  dynamically  increasing  load  resistance  on  normalized  power 
output  relative  to  normalized  load  resistance. 


4.5.  Dynamic  Load  resistance  on  internal  resistance  of  the  generator 

The  influence  of  an  increasing  electrical  load  is  investigated  in 
order  to  verify  the  postulation  that  electrical  load  matching  yields 


Fig.  15.  R[  for  varying  thermal  input  conditions. 


peak  thermoelectric  power  under  dynamic  conditions.  Indeed,  the 
maximum  point  on  the  curve  generated  by  Eq.  (7)  occurs  when 
RL  =  R int-  To  this  end,  from  circuit  theory,  the  internal  resistance 
of  the  generator  is  calculated  by  measuring  the  open  circuit  voltage 
[20,23]: 

Rm  =  ^-Rl.  (16) 

The  validity  of  electric  load  matching  for  peak  power  under  an 
increasing  load  is  examined  under  conditions  in  which  the  thermal 
input  varies  with  varying  hot  fluid  flow,  cold  fluid  flow  and  tem¬ 
perature.  The  results  are  illustrated  in  Fig.  14.  The  flow  rates  are  ex¬ 
pressed  in  L/min  and  the  temperatures  in  °C.  The  electrical  load 
resistance  is  increased  at  a  rate  of  5.8  ±  2.2  Q/s. 

Fig.  14  shows  that,  for  all  thermal  conditions  tested,  electrical 
load  matching  yields  a  minimum  of  97.6%  peak  power.  The  near 
peak  performance  is  attributed  to  the  fact  that  the  measured  RLy0pt 
is  less  than  the  internal  resistance  of  the  generator.  For  this  reason, 
the  electrical  load  matching  implies  a  load  resistance  that  is  neces¬ 
sarily  greater  than  optimal  which,  as  illustrated  in  Fig.  1 1,  is  a  more 
forgiving  error  than  an  electrical  load  that  is  less  than  optimal;  this 
is  due  to  a  relatively  small  power  rate  of  change  with  respect  to 
electrical  load.  Indeed,  it  is  shown  in  Fig.  15  that  the  optimal  value 
of  R[  is  strictly  less  than  one  ranging  from  0.76  to  0.92. 

The  results  presented  in  Fig.  15  are  consistent  with  the  observa¬ 
tions  made  in  [13]  on  a  single  thermoelectric  module  in  which  the 
optimal  electrical  load  produced  from  Eq.  (16)  was  strictly  greater 
than  the  measured  optimal  load. 


5.  Conclusion 

A  forty  module  thermoelectric  liquid-to-liquid  generator  is  used 
to  measure  thermoelectric  power  production  under  an  increasing 
electrical  load  resistance.  Varying  thermal  conditions  show  that 
power  output  behaves  exponentially  relative  to  temperature  and 
as  a  root  function  relative  to  fluid  flow.  The  thermal  input  condi¬ 
tions  are  quantified  in  terms  of  a  non-dimensional  term  that  incor¬ 
porates  flow  rates  and  temperatures.  A  correlation  is  developed 
predicting  peak  power  output  to  within  a  10%  error  margin.  The 
optimal  electrical  load  resistance  remained  within  15%  of  15.8  Cl 
for  all  thermal  input  conditions  tested  yielding  a  minimum  of 
99.0%  of  the  measurable  peak  power  and  making  an  MPPT  unnec¬ 
essary  for  many  applications.  Furthermore,  it  is  shown  that  ther¬ 
moelectric  power  production  is  less  sensitive  to  electrical  load 
error  when  the  load  is  greater  than  optimal  rather  than  less  than 
optimal.  Since  it  is  also  shown  that  electrical  load  matching  consis¬ 
tently  generates  electrical  loads  that  are  greater  than  optimal,  load 
matching  is  deemed  in  a  zone  less  sensitive  to  error  yielding  97%  or 
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more  of  the  measurable  peak  power  for  all  thermal  input  condi¬ 
tions  tested. 
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